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We present the results of a set of Monte Carlo simulations of Dynamically Triangulated Random Surfaces 
embedded in three dimensions with an extrinsic curvature dependent action. We analyze several observables 
in the crossover regime and discuss whether or not our observations are indicative of the presence of a phase 
transition. 
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1. Introduction 

In this work, we investigate a theory of fluid, 
fluctuating random surfaces embedded in three 
dimensions. Theories of fluctuating surfaces 
(string theories) have been conjectured to de- 
scribe a wide variety of physical phenomena and 
models, including the strong interaction at large 
distances, the 3d Ising model, and unified the- 
ories incorporating gravity. Lipid bilayers and 
microemulsions are also examples of fluctuating 
surfaces [0. Such biological and chemical mem- 
branes exhibit self-avoidance, which we do not 
take into account in our simulations. Models of 
fluctuating random surfaces can in fact be solved 
exactly when the surfaces are embedded in dimen- 
sions D < 1; these solutions break down though 
when continued to the more physical regime D > 
1. This may be related to the observation that 
simulations of fluctuating surfaces in Z? > 1 using 
a discretization of the standard Polyakov string 
action are dominated by crumpled, spiky config- 
urations. 

Our lattice model is constructed by triangulat- 
ing each surface. Each node of the triangulation 
is embedded in by the functions Xf ; i labels 
the ith node and /i runs from 1 to 3. The tri- 
angulation is characterized by the adjacency ma- 



trix Cij, whose elements equal 1 if i and j label 
neighboring nodes of the triangulation and van- 
ish otherwise. The triangles are assumed to be 
equilateral (as measured by the intrinsic metric 
of each surface); the connectivity at each node 
determines the intrinsic curvature. We also as- 
sociate a normal vector {n^)f. with each triangle 
(indices with hats label the triangles). We shall 
study the theory defined by the action 
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For A = this is simply a discretization of the 
Polyakov string action. The final term, which de- 
pends on the discretized extrinsic curvature, in- 
troduces a ferromagnetic interaction between sur- 
face normals, which one might hope would cause 
smoother surfaces to dominate the partition func- 
tion. We would like to know if there is a smooth 
phase and a phase transition (at some finite Ac) 
between this phase and the crumpled phase ob- 
served at A = 0. If this were so, an interest- 
ing continuum limit of this lattice model could 
perhaps be constructed at this phase transition 
point, yielding a new continuum string theory. 

2. The Simulation 

We have considered triangulations with the 
topology of the torus, to minimize finite size ef- 
fects. The above action was used, with the BRST 



2 



invariant measure utilized also by Baillie, John- 
ston and Williams so that 

Z= j\{dXt\{qfeM~SG-\SE).{2) 

GGT(l) i 

where D = 3, qi is the connectivity of the ith 
vertex and T(l) refers to the set of triangulations 
of genus 1. We used the standard Metropolis al- 
gorithm to update our configurations. To sweep 
through the space of triangulations we performed 
flips on randomly chosen links. Flips were au- 
tomatically rejected if they yielded a degenerate 
triangulation. After a set of 3M flips was per- 
formed, 3M randomly selected embedding coor- 
dinates were updated via random shifts from a 
flat distribution. Most of the Monte Carlo simu- 
lations were performed on HP-9000 (720 and 750 
series) workstations; we also collected some data 
by simulating lattices on each of the 32 nodes of 
a CM-5. 

We ran on lattices ranging in size from iV = 36 
to 576 (N signifles the number of vertices) with 
4 to 7 different values of A for each N. Most of 
the data was this data was taken in the region 
A G (1.325, 1.475). For smaU TV, the runs con- 
sisted of 3 X 10^ sweeps, while we performed longer 
runs (of up to 27 x 10^ sweeps for TV = 576) for 
larger lattices, because the auto-correlation times 
for our simulations were very large. (The corre- 
lation time for the radius of gyration was greater 
than 10^ sweeps for N = 576!) To determine our 
observables as a function of A we used a histogram 
reconstruction procedure. We patched different 
histograms by weighting them with the associ- 
ated statistical indetermination (which was esti- 
mated by a jack-knife binned procedure). Various 
consistency checks indicate that this procedure is 
very reliable. 

3. Observables 

We measured the edge action Se and the as- 
sociated specific heat C(A) = ^(< S% > — < 
Se >^)- In Fig- 1 we plot the specific heat curve 
(constructed via the histogram procedure) and we 
tabulate its maximum and peak position for var- 
ious lattice sizes in Table 1. 



Figure 1. C(A). Dotted fines: N=144. Dashed 
lines: N=288. Solid lines: N=576. 
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36 


3.484(8) 


1.425(35) 


72 


4.571(15) 


1.410(15) 


144 


5.37(14) 


1.395(30) 


288 


5.55(7) 


1.410(25) 


576 


5.81(17) 


1.425(30) 



Table 1 

The maximum of the specific heat and its posi- 
tion, with errors, for different lattice sizes. 



We see that the specific heat peak grows vig- 
orously with N for small lattices, but that this 
growth quickly levels off for larger N. These ob- 
servations agree fairly well with previous work |^ 

. For the larger lattices it appears that the peak 
position shifts very slowly towards higher values 
of A, though this increase is not statistically sig- 
nificant. The shape of the peak does not change 
dramatically with N] it narrows perhaps a bit 
between TV = 144 and TV = 576. 

We measured the magnitude of the extrinsic 
Gaussian curvature, J \ K \ y/\ h \ {h is the in- 
duced metric and K is the determinant of the 
extrinsic curvature matrix), given by 
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Here (f)-[ denotes the angle subtended by the jth 
triangle at the ith vertex. This quantity, plotted 
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in Fig. 2, measures the magnitude of the deficit 
angle in the embedding space averaged over all 
vertices. Note that the mean Gaussian curvature 



Figure 2. | /C |, plotted as in Figure 1. 



decreases rapidly in the neighborhood of A = 1.4, 
indicating that a sharp crossover is occuring in 
this system. From this plot we can see that finite 
size effects increase with A; they do not appear 
to peak in the region about A = 1.4 as one might 
expect for a typical phase transition. We also 
measured various other observables which char- 
acterize both the intrinsic and extrinsic geometry 
of these surfaces (and the correlation between in- 
trinsic and extrinsic geometry). These measure- 
ments are discussed in a longer write-up of this 
work 1^. They all exhibit sharp crossover be- 
havior in the region near A = 1.4. We found 
that the auto-correlation times of these observ- 
ables grew rapidly as A increased, but we did not 
note any maximum in these times in the region 
about A = 1.4. 

4. Interpretation 

This model of crumpled surfaces appears to 
exhibit sharp crossover behavior in the region 
around A = 1.4. The sharp change in the mag- 
nitude of the Gaussian curvature, the radius of 
gyration and other observables indicates that the 
normals acquire long-range correlations, up to 
the size of the systems we examine. The zero 



string tension measurement of also shows that 
the disordered regime differs from the regime in 
which the surfaces are ordered (up to scale of 
the lattices that are simulated) by only a small 
shift in A. This evidence might indicate the pres- 
ence of a phase transition at this point. Since 
the peak growth rapidly diminishes for large iV, 
such a phase transition would likely be higher 
than second order. Still, the apparent absence of 
diverging correlation times and increasing finite 
size effects in the peak region leads us to ques- 
tion whether we are actually observing a typical 
phase transition. 

There are indeed other possible interpretations 
of our data. Note that the surfaces which we sim- 
ulate are quite small. For instance, if the surfaces 
in our simulations had an intrinsic dimension of 
2.87 (characteristic oi D = gravity), they would 
have roughly a linear size of fewer than 10 lattice 
spacings. 

Perhaps the simplest alternative explanation 
for the presence of this peak is suggested by the 
arguments of Kroll and Gompper |9| . They argue 
that the peak occurs when the persistence length 
of the system approaches the size of the lattice 
(Cj9 ~ ^^); d denotes the intrinsic Hausdorff di- 
mension. Fluctuations on a larger scale become 
more important. When this scale is greater than 
the lattice size these fluctuations are suppressed. 
Thus one might surmise that the specific heat 
will drop for large A. The one-loop renormal- 
ization group calculation jl^ predicts that the 
persistence length grows as ~ exp(CA); C is 
inversely proportional to the leading coefficient 
of the beta function. We would expect that the 
peak position should shift to the right with in- 
creasing N in this scenario as 



\eak{N') — Xpeak{N) 



ln(f) 
dC ' 



(4) 



This reasoning also indicates that the peak should 
widen as the lattice size increases; we do not ob- 
serve this at all. 

An alternative scenario, which builds on the 
ideas in the above paragraph, is suggested by the 
tantalizing similarities between the results of our 
fluid surface simulations and what has been ob- 
served for the d = 4 SU{2) Lattice Gauge Theory 
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||n| and for the d = 2 0(3) model. 

The 0(3) model, which is thought to be asymp- 
totically free, exhibits a specific heat peak near 
/? = 1.4 (first measured via Monte Carlo simu- 
lations by Colotjll). The origin of this peak 
is understood |13|; it is due to the fluctuations 
of the sigma particle, a low-mass bound state 
of the massless 0(3) pions. The sigma induces 
short-range order and contributes to the specific 
heat as a degree of freedom only at high temper- 
atures (when the correlation length in the sys- 
tem becomes smaller than its inverse mass). The 
peak thus occurs at the beginning of the crossover 
regime, when the correlation length is several lat- 
tice spacings. 

According to the low temperature expansion, 
the correlation length grows as ^ ~ exp(27r /?)//?. 
Thus one would expect a fairly rapid crossover in 
the 0(3) model; the correlation length should in- 
crease by roughly a factor of 9 when fi is shifted 
by about 0.35. Such a crossover is indeed ob- 
served, though it is not so apparent that it is as 
dramatic as the crossover behaviour observed for 
fluid surfaces. 

Recent simulations of the 0(3) model in- 
dicate that the speciflc heat peak grows signifl- 
cantly when the lattice size L is increased from 5 
to 15, and that virtually no growth in peak height 
is evident as L is increased further up to 100. Fur- 
thermore, the peak position shifts to the right as 
L grows and then appears to stabilize for large 
L. This is more or less what we observe in our 
simulations of fluid surfaces, on lattices of small 
size. We point out these similarities largely to 
emphasize that there does exist an asymptotically 
free theory (with low mass excitations) which ex- 
hibits crossover behavior qualitatively similar to 
that observed in our simulations. The analogy 
is perhaps deeper, though, since the fluid surface 
action (with extrinsic curvature) in certain guises 
looks like a sigma model action. It would not 
therefore be so surprising from this point of view 
to flnd a sigma particle in these theories, perhaps 
associated with {ti? — 1), in which fi denotes the 
unit normal to our surfaces. 
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